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Many fractional quantum Hall wave functions are known to be unique highest-density zero modes 
of certain “pseudopotential” Hamiltonians. While a systematic method to construct such parent 
Hamiltonians has been available for the infinite plane and sphere geometries, the generalization to 
manifolds where relative angular momentum is not an exact quantum number, i.e. the cylinder 
or torus, has remained an open problem. This is particularly true for non-Abelian states, such as 
the Read-Rezayi series (in particular, the Moore-Read and Read-Rezayi Z 3 states) and more exotic 
non-unitary (Haldane-Rezayi, Gaffnian) or irrational (Haffnian) states, whose parent Hamiltoni¬ 
ans involve complicated many-body interactions. Here we develop a universal geometric approach 
for constructing pseudopotential Hamiltonians that is applicable to all geometries. Our method 
straightforwardly generalizes to the multicomponent SU(n) cases with a combination of spin or 
pseudospin (layer, subband, valley) degrees of freedom. We demonstrate the utility of our approach 
through several examples, some of which involve non-Abelian multicomponent states whose par¬ 
ent Hamiltonians were previously unknown, and verify the results by numerically computing their 
entanglement properties. 

PACS numbers: 73.43.-f 


I. INTRODUCTION 

Among the most striking emergent phenomena in con¬ 
densed matter are the incompressible quantum fluids in 
the regime of the fractional quantum Hall effectP^. A key 
theoretical insight to understanding the many-body na¬ 
ture of such phases of matter was provided by Laugh- 
lin’s wave functioiP. Shortly thereafter, Haldan^ re¬ 
alized that the = 1/3 Laughlin state also occurs as 
an exact zero-energy ground state of a certain positive 
semi-definite short-range parent Hamiltonian which an¬ 
nihilates any electron pair in a relative angular momen¬ 
tum I > 1 state, but which assigns an energy penalty for 
each state with / = 1. 

Haldane’s construction of the parent Hamiltonian for 
v = 1/3 is just one aspect of a more general “pseu¬ 
dopotential” formalism that applies to quantum Hall 
systems in the infinite plane or on the surface of a 
spher^. In those cases, the system is invariant under 
rotations around at least a single axis, and by virtue of 
the Wigner-Eckart theorem, any long range interaction 
(such as Coulomb interactions projected to Landau level) 
decomposes into a discrete sum of components VJ. The 
different components are quantized according to the rela¬ 
tive angular momentum /, which is odd for spin-polarized 
fermions and even for spin-polarized bosons. The unique 
zero mode of the Vi pseudopotential at = 1/3 is pre¬ 
cisely the Laughlin state. Furthermore, it is the densest 
such mode because all the other states, at u = 1/3 or 
any filling C > n, are separated by a finite excitation 
gap as indicated by overwhelming numerical evidence. 
In numerical simulations, it is furthermore possible to 


selectively turn on the magnitude of longer range pseu¬ 
dopotentials (V 3 ,V 5 , etc.), and verify that the Laughlin 
state adiabatically evolves to the exact ground state of 
the Coulomb interactioiP. The corrections induced to 
the Laughlin state in this way are notably small (below 
1% in finite systems containing about 10 particles), and 
the gap is maintained during the proces^. These findings 
constitute an important support of Laughlin’s theory. 

Remarkably, the existence of pseudopotential Hamil¬ 
tonians is not limited to the Laughlin states. Gener- 
ically, more complicated parent Hamiltonians arise in 
the case of quantum trial states with non-Abelian quasi¬ 
particles, which are one important route towards topo¬ 
logical quantum computatiorl^. For example, the cel¬ 
ebrated Moore-Read “Pfaffian” stat^, believed to de¬ 
scribe the quantum Hall plateau at n = 5/2 filling frac- 
tiorP, was shown to possess a parent Hamiltonian which 
is the shortest-r ang e repulsive potential acting on 3 par¬ 
ticles at a time^Hi^. This state is a member of a family of 
states - the parafermion Read-Rezayi (RR) sequence - 
where the parent Hamiltonian of the kth member is the 
shortest-range /c-body pseudopotentiaP^. Similar par¬ 
ent Hamiltonians have been formulated for non-Abelian 
S = k/2 chiral spin liquid stateMore recent analy¬ 
sis has allowed to resolve the structure of quantum Hall 
trial states interpreted as null spaces of pseudopotential 
HamiltonianJ^^l®. 

The knowledge of the parent Hamiltonian class is cru¬ 
cial for the complete characterization of a quantum Hall 
state. In certain cases, when the state can be represented 
as a correlator of a conformal field theory (CFT), the 
charged (quasihole) excitations can be constructed using 
the tools of the CFTP. (Note that such approaches be- 
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come rather cumbersome for quasielectron excitationJI^, 
and even more so in the torus geometr}!^.) The tools 
of CFT relating wave functions to conformal blocks cer¬ 
tainly become less useful when information about neu¬ 
tral excitations is needed. The knowledge of the par¬ 
ent Hamiltonian is thus indispensible, e.g. for estimating 
the neutral excitation gap of the system. In some cases, 
the neutral gap c an be computed by the single-mode ap- 
proximatiorPi^, which is microscopically accurate for 
Abelian stateP^, but requires non-trivial generalizations 
for the non-Abelian stateP^. Similarly, entanglement 
spectra may help to characterize the elementary exci¬ 
tations solely deduced from the ground state wave func- 
tioiP^Hm^ but only unfold their full strength as a com¬ 
plementary tool to the spectral analysis of the associated 
parent Hamiltonian. While a CFT underlies the struc¬ 
ture of entanglement, conformal blocks, and clustering 
properties within most quantum Hall states of interest!^, 
it is desirable to have a general framework that comple¬ 
ments it with a pseudopotential parent Hamiltonian. 

An appealing “added value” of pseudopotentials as 
building blocks of parent Hamiltonians is the possibil¬ 
ity of discovering new states by varying the considered 
pseudopotential terms and searching for new zero-mode 
ground states. This strategy is exemplified by the spin- 
singlet Haldane-Rezayi stat^^, initially discovered as a 
zero mode of the “hollow core” interaction between spin¬ 
ful fermions at half filling. Another example is the 3-body 
interaction of a slightly longer range than the one that 
gives rise to the Moore-Read state. It was founcP^ that 
such an interaction has the densest zero-energy ground 
state at filling factor u = 2/5, subsequently named the 
“Gaffnian’®. These examples illustrate that a system¬ 
atic description of an operator space of parent Hamil¬ 
tonians holds the promise of the discovery of previously 
unknown quantum Hall trial states. 

All appreciable aspects of parent Hamiltonians dis¬ 
cussed so far are independent of the geometry of the man¬ 
ifold in which the quantum Hall state is embedded. For 
several purposes, however, knowing the parent Hamilto¬ 
nian on geometries without rotation symmetry, such as 
the torus or cylinder, is particularly desirable. For exam¬ 
ple, quantum Hall states only exhibit topological ground 
state degeneracy on higher genus manifolds such as the 
toru^^. Accessing the set of topologically degenerate 
ground states further allows to extract the modular 5, T 
matrices which e ncode all topological information about 
the quasiparticle^^^^. Furthermore, parent Hamiltoni¬ 
ans on a cylinder or torus can be used to derive solv¬ 
able models for quantum Hall states when one of the 
spatial dimensions becomes comparable to the magnetic 
lengtlpUM]^ It has been shown that such models can be 
used to construct “matrix-product state” representations 
for quantum Hall states, and in some cases can be used 
to study the physics of “non-unitary” states and classify 
their gapless excitation^^^^^. 

A systematic construction of many-body parent Hamil¬ 
tonians for quantum Hall states was first undertaken by 


Si mon et al for the infinite plane or sphere geome- 
tr^ F^F^ i This approach relies on the relative angular mo¬ 
mentum which in this case is an exact quantum number. 
As such, it cannot be applied to the cylinder, torus, or 
any quantum Hall lattice model such as fractional Chern 
insulators, for which the analogue of pseudopotentials has 
been developed recent l}®HiII^ 

Ref. [44] has introduced the closed-form expressions for 
all two-body Haldane pseudopotentials on the torus and 
cylinder. In this work, inspired by Refs. |42]-[44] as the 
starting point of our analysis, we provide a complete 
framework for constructing general quantum Hall par¬ 
ent Hamiltonians involving Wbody pseudopotentials, for 
fermions as well as bosons, in cylindric and toroidal ge¬ 
ometries. This advance proves particularly important for 
the non-Abelian states, most of which necessitate many- 
body pseudopotentials in their parent Hamiltonian class. 
From the construction scheme laid out in this work, all 
topological properties of the non-Abelian states such as 
their modular matrices and topological ground state de¬ 
generacy can now be conveniently studied from their as¬ 
sociated toroidal parent Hamiltonian. Complementing 
previous results for the sphere and infinite plane, our for¬ 
malism furthermore directly generalizes to multicompo¬ 
nent systems with an arbitrary number of “spin types” or 
“colors”. Therefore, our construction of many-body clus¬ 
tered Hamiltonians not only applies to arbitrary geome¬ 
tries, but also crucially simplifies previous approaches. 
We illustrate this by numerical examples, including non- 
Abelian multicomponent states whose parent Hamiltoni¬ 
ans were previously unknown. 

The article is organized as follows. In Sec. |TI| we start 
with a brief overview of Haldane’s pseudopotential for¬ 
malism from the viewpoint of clustering conditions, and 
define our notation of generalized many-body multicom¬ 
ponent pseudopotentials. The clustering conditions, as 
well as the polynomial constraints following from them, 
form the core of our systematic geometric construction 
of many-body parent Hamiltonians described in Sec. |ln| 
The main result of that section is the appropriately cho¬ 
sen integral measure which allows us to generate all pseu¬ 
dopotentials through a direct Gram-Schmidt orthogonal- 
ization. The extension to many-body parent Hamiltoni¬ 
ans for spinful states is described in Sec.[lV| Technical de¬ 
tails of the construction are delegated to the appendices. 
In Section |V| explicit examples of parent Hamiltonian 
studies are worked out in detail. We discuss spin polar¬ 
ized states such as the Gaffnian, the 1/q Pfaffian and the 
Haffnian states, as well as spinful states such as the spin- 
singlet Gaffnian, the NASS and the Halperin-permanent 
states. Apart from deriving the second-quantized par¬ 
ent Hamiltonians, we also provide extensive numerical 
checks of our construction using exact diagonalization 
and analysis of entanglement spectra. In Sec. jVI] we con¬ 
clude that our geometric construction of parent Hamilto¬ 
nians promises ubiquitous use in the analysis of fractional 
quantum Hall states and outline a few immediate future 
directions. 
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II. HALDANE PSEUDOPOTENTIALS 


Within a given Landau level, the kinetic energy term 
in the quantum Hall (QH) Hamiltonian is “quenched”, 
i.e. is effectively a constant. Hence the remaining ef¬ 
fective Hamiltonian only depends on the interaction be¬ 
tween particles (e.g., the Coulomb potential) projected to 
the given Landau leveP. In the infinite plane, the lowest 
Landau level (LLL) projection amounts to evaluating the 
matrix elements of the Coulomb interaction between the 
single-particle states of the form 

(pm{z) = z"^ exp{-\z\'^ (1) 

where ^ is a complex parametrization of 2D electron coor¬ 
dinates and Ib = \/h/eB the magnetic length (Fig. [^). 
For simplicity, we consider a QH system in the back¬ 
ground of a fixed (isotropic) metric^, which allows us 
to write ^ = x iy. The states in Eq. are mutually 
orthonormal, and span the basis of the LLL. There are 
of these states, which is also the number of mag¬ 
netic fiux quanta through the system. The above will 
be assumed throughout this paper, which is appropriate 
in strong magnetic fields when the particle-hole excita¬ 
tions to other Landau levels are suppressed by the large 
cyclotron energy gap hwc. 
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FIG. 1. (Color online) (a) Landau levels and the single¬ 
particle states on the disk. The lowest LL is separated from 
the next lowest level by the cyclotron gap hojc- The wave 
functions of the LLL states are (j)rn^ labelled by an integer 
m = 0,1, 2..., and Gaussian-localized along the ring of ra¬ 
dius ~ (b) Any two-particle state separates into a 

product of two wave functions, one depending on the center- 
of-mass and the other depending on the relative coordinate 
z = zi — Z 2 . The relative wave function 'i^rn has identical form 
to the single-particle wave function, with two important differ¬ 
ences: its effective magnetic length is rescaled {is ^ 
and the integer value of m is constrained by particle statistics 
(assuming the spin is fully polarized): for fermions, the m’s 
entering the relative wave function are odd (red dashed lines). 


(Fig#). As usual, the solution of the two-body problem 
proceeds by transforming from coordinates Zi , Z 2 into the 
center of mass (COM) Z = [ziZ 2)/2 and the relative 
coordinate z = zi — Z 2 frame. In the new coordinates, 
the two-particle wave function decouples. As we are in¬ 
terested in translationally-invariant problems, only the 
relative wave function (which depends on z) will play a 
fundamental role in the following analysis. For any two 
particles, the relative wave function turns out to have an 
identical form to the single-particle wave function 0 

ipm{z = zi- Z2) cx 2:™ exp(-|2;|^/8£|), (2) 

up to the rescaling of the magnetic length is- An im¬ 
portant difference between Eqs. 0 and 0 is the new 
meaning of m: since ^ now represents the relative sepa¬ 
ration between two particles, m in Eq. ([^ is related to 
particle statistics, and therefore encodes the clustering 
properties. Eor spinless electrons, m in Eq. 0 is only al¬ 
lowed to take odd integer values since the wave function 
must be antisymmetric with respect to ^ — z, while 
for spinless bosons m can be only be an even integer. 
Einally, m is also the eigenvalue of the relative angular 
momentum for two particles (h = 1), as we can directly 
confirm from 


i 

After this two-particle analysis (summarized in Eig.[^, 
we are in position to introduce the notion of clustering 
properties for A-particle states. Let us pick a pair of 
coordinates Zi and Zj of indistiguishable particles in a 
many-particle wave function. We say that these particles 
are in a state which obeys the clustering property 
with the power m if vanishes as a polynomial of total 
power m as the coordinates of the two particles approach 
each other: 


■tpilizi Zj) ~ {zi - Zj)'^. (4) 

Similarly as before, we can relate the exponent m to the 
angular momentum if the latter is a conserved quan¬ 
tity. Clustering conditions like this directly generalize to 
cases where more than 2 particles approach each other, 
with the polynomial decay also specified by a power. Eor 
example, we say that an A-tuple of particles 2 : 1 , 2 ( 2 ,... 2 :at 
is in a state with total relative angular momentum m 
if -0^ vanishes as a polynomial of total degree m as the 
coordinates of N particles approach each other. Eixing 
an arbitrarily chosen reference particle of the A-tuple, 
i.e. zi =: Zr^ we have: 


Restricting to a single LL, a large class of QH states can 
be classified by their elustering properties. (See Refs. [49] 
and EOl for classification schemes based on clustering.) 
These are a set of rules which describe how the wave func¬ 
tion vanishes as particles are brought together in space. 
To define the clustering rules, it is essential to first con¬ 
sider the problem of two particles restricted to the LLL 


N 

i^^{Z2,...,ZN Zr) ~ 0(2^ - Zj)^\ (5) 

i=2 

with ^ as all remaining N — 1 particles ap¬ 

proach the reference particle Zr. If the system is rota- 
tionally invariant about at least a single axis (such as 
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for a disk or a sphere) it directly follows that the state 
in Eq. is also an eigenstate of the corresponding Lz 
relative angular momentum operator with eigenvalue m. 

The simplest illustration of the clustering condition is 
the fully filled Landau level. The wave function for such 
a state is the single Slater determinant of states (j)m{zi) 
in Eq. Due to the Vandermonde identity, this wave 
function can be expressed as 

~ ^ 3 ) exp(- y] \zk\^/Al%). (6) 

i<j k 

We see that when any pair of particles Zi and Zj is iso¬ 
lated, the relevant part of the wave function is {zi — Zj). 
Therefore, the wave function of the filled Landau level 
vanishes with the exponent m = 1 as particles are 
brought together. This is the minimal clustering con¬ 
straint that any spinless fermionic wave function must 
satisfy. (As we will see below, interesting many-body 
physics results from stronger clustering conditions on the 
wave function.) 

When properly orthogonalized, the states form 

an orthonormal basis in the space of magnetic translation 
invariant QH states. Th^ allow to define the V-body 
Haldane pseudopotential^ (PPs) 

= I'lpm) {^m\ , (7) 

which obey the null space condition = 0 for 

m' 7 ^ m. Since they are positive-definite, the PPs 
give energy penalties to V-body states with total relative 
angular momentum m. With a given many-body wave 
function, the Hamiltonian representation of will in¬ 
volve the sum over all V-tuple subsets of particles. Con¬ 
crete examples of this will be given in Secs. |In]and|Vl 

Eor a given filling fraction, it can occur that a certain 
QH state is the unique and densest ground state lying 
in the null space of a certain linear combination of PPs, 
i.e. it is annihilated by a certain number of PPs. (The 
requirement of being the densest state is necessary to ren¬ 
der the finding non-trivial, because it is simple in prin¬ 
ciple to construct additional zero modes of a given PP 
Hamiltonian by increasing the magnetic flux, i.e. by nu¬ 
cleating quasihole excitations.) The most elementary ex¬ 
amples are the Laughlin states at 1/m filling, which lie in 
the null space of U‘^, for all m' < m. As it also represents 
the densest configuration that is annihilated by the PP, 
the Laughlin state emerges as the unique ground state 
of a Hamiltonian at filling 1/m, H = Yjm'<m^rn'U^,, 
where the coefficients are arbitrary as long as Cm' > 0. 
That is, the fermionic 1/3 Laughlin state is the unique 
ground state of , while the fermionic 1/5 state is the 
unique groudstate of any linear combination of Ui and 
t /3 with positive weights. Note that the PPs of even 
m' are precluded by fermionic antisymmetry. Eor short- 
ranged two-body interactions on the disk where the de¬ 
gree of the clustering polynomial m coincide with the 
exact relative angular momentum /, the traditional no¬ 
tation by Haldan^ relates to ours via Vi = U^. As we 


elaborate below, many more exotic states can be realized 
as the highest density ground states of combinations of 
PPs involving N >2 bodies. 

In view of other geometries than disk or sphere, one 
question immediately arises: Is there any hope of defin¬ 
ing PPs in the ahsenee of eontinous rotation symmetry 
and hence no exact relative angular momentum quan¬ 
tum number? This question is natural because one of 
the popular choices for the gauge of the magnetic field 
- the Landau gauge - is only compatible with periodic 
boundary conditions along one or both directions in the 
plane, which breaks continuous rotational symmetry. 

The answer to the above question is affirmative^ which 
is the central message of this paper. This is mainly be¬ 
cause the clustering conditions, i.e. the polynomial ex¬ 
ponent m, are more fundamental than their interpreta¬ 
tion as the relative angular momentum. The clustering 
conditions are short-distance properties: their support is 
the area associated with the fundamental droplet of N 
particles ^ N x 27r^^. Assuming a one-to-one correspon¬ 
dence between a clustering power and a pseudopotential, 
the PPs should be independent of the specific geome¬ 
try as long as they act in a manifold that is homoge¬ 
neous and much larger than the fundamental A'-particle 
droplet. This viewpoint has been confirmed by the ex¬ 
plicit constructions of quantum Hall trial wave functions. 
Eor instance, the successful generalization of the Laugh¬ 
lin wave function to the torus, given in the classic pa¬ 
per by Haldane and RezayJ^, has demonstrated that its 
short-distance properties are identical to its original ver¬ 
sion defined on the disk. The states in both geometries 
are uniquely characterized by their clustering properties 
and are locally indistinguishable; their main difference 
lies in the global properties, e.g., the fact that the torus 
Laughlin state is m-fold degenerate due to its invariance 
under the center of mass (COM) translation. This de¬ 
generacy is of intrinsically topological origiiP^. At fill¬ 
ing u = pjq^ where p^q are relatively co-prime integers, 
the quantum Hall state on a torus is invariant under the 
translation that moves every particle by q orbitals (see 
Fig-i; below). This symmetry guarantees an exact m- 
fold degeneracy of the Laughlin u = 1/m state on the 
toru^. Additional topological degene racie s can arise for 
more complicated (non-Abelian) state^^^. 

In the following, we assume there exists, in general, a 
well-defined deformation of the null space specified by a 
planar QH trial wave function to the multi-dimensional 
null space specified by the associated set of topologi¬ 
cally degenerate QH ground states on the torus. What 
we are then interested in is to find a suitable deforma¬ 
tion of the planar Laplacian, whose bilinear form is the 
known planar parent Hamiltonian composed of the spher¬ 
ical PPs, to the toroidal Laplacian, whose bilinear form 
is the toroidal Hamiltonian. In the following sections, we 
solve this problem by what we refer to as the “geometric 
construction” of pseudopotentials. 
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III. GEOMETRIC CONSTRUCTION OF 
PSEUDOPOTENTIALS: SPINLESS CASE 


We describe the construction of generic QH pseudopo¬ 
tential Hamiltonians with a perpendicular magnetic field 
applied to the 2D electron gas. We assume the field to 
be sufficiently strong such that the spin is fully polarized, 
and that there are no further internal degrees of freedom 
for the particles. We first introduce a suitable single¬ 
particle basis for generic A-body interactions that obeys 
the magnetic translation symmetry and conserves the 
center-of-mass (COM) momentum. Next, we show how 
the explicit functional form of Haldane PPs in Sec.[ll|can 
be easily obtained from geometric principles and symme¬ 
try. We constrain ourselves to single-component PPs in 
this section, and generalize our construction to multicom¬ 


ponent (spinful) PPs in Sec. IV 


A. Basis choice 

A well-developed p seud opotential formalism is avail¬ 
able in the literatur^^^^ for the infinite plane or the 
sphere, where the 2 ;-component of angular momentum 
is conserved. A different approach to pseudopotential 
Hamiltonian construction, however, is needed when the 
system is no longer invariant under continuous rotations 
around the z-axis. This can occur when periodic bound¬ 
ary conditions are imposed [Fig|2], either along one di¬ 
rection (cylinder geometry) or both directions (torus). 
Under periodic boundary conditions (PBCs) in one direc- 


n = 2'KlBlLy 






ny—\{x I i B — 


^r-(r) = Efc^^n+fciv,(r) 


FIG. 2. (Color online) Boundary conditions considered in 
this paper, compatible with Landau gauge: cylinder (left) and 
torus (right). The associated single-particle states '0n(r) are 
arranged along a ID chain, with the parameter n — 
controlling the distance between nearest neighbor orbitals. 


tion (say the single-particle Hilbert space is spanned 
by the Landau gauge basis wave function labelled by n: 

V’n(a^,2/) = (F|c], |0) oc , (8) 

where Ib = is the magnetic length and cj^ is 

the second-quantized operator that creates a particle in 
the state |n). The parameter k — 2'KiBlLy sets the ef¬ 
fective separation between the one-body states in the x- 
direction, as each one-body state is a Gaussian packet ap¬ 
proximately localized around K,n in x-direction (Fig. [^. 

Due to this simple one-parameter labeling of the one- 
body states, an V-body interaction matrix element is la¬ 


beled by 2N indices. Additional constraints on its func¬ 
tional form arise when we consider interactions projected 
to a given Landau level (Fig. |^. Due to the magnetic 
translation invariance, the interactions projected to a 
Landau level only allow for scattering that conserves to¬ 
tal momentum, or equivalently leaves the COM of the 
particles fixed. For example, the process in Fig. IT is 
allowed, but the one in Fig.[^ is forbidden. This special 
structure in the interaction Hamiltonian gives rise to the 
symmetry under many-body translations that shifts ev¬ 
ery particle by q orbitals at filling u = pfq (Fig. §)• This 
is the symmetry that underlies the topological ground 
state degeneracy. 


(a) 



(b) > 



r 1 Q 



FIG. 3. (Golor online) Scattering schemes and symmetries 


in a Landau level. Preservation of GOM allows a cluster of 
particles to scatter according to (a), but forbids the scattering 
according to (b). (c) The many-particle translation operator 


T acts on a given configuration by moving every particle q 
orbitals to the right (in this example, q = 2). Many-body 
states at filling v — pjq are invariant under this symmetry, 
which gives rise to the topological degeneracy equal to q. 


The matrix elements of any interaction U(ri,..., rjv) 
projected to a Landau level are given by 

= j dri...drN 

r(ri,...,riv), (9) 

which corresponds to the second-quantized Hamiltonian 


a n 

( 10 ) 

Here fij = rij—n/N denotes the orbital of particle j with 
respect to the COM, n = ni + ... + nvr, and Poc is a poly¬ 
nomial in N variables Bfij. The index a, as will be clear 
from the explicit construction of Pa below, specifies the 
degree of the polynomial and, for a multicomponent state 
discussed in Section [rV| its spin sector. The polynomial 
Pa can be chosen to be real, as will be evident from its 
geometric construction to follow. From now on, n will 
refer to the COM, and not the index of a single particle 
previously appearing in Eq. c\^. is the same operator 
as in Eq. creating the jth particle in state \nj). The 
Hamiltonian thus consists of a product sum over all po¬ 
sitions of the COM n, as well as all polynomial degrees 
a. 
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Eqs. 1^ and represent a general translationally in¬ 
variant Hamiltonian projected to a Landau level. This 
Hamiltonian has a rather special form: it decomposes 
into a linear combination of positive-definite operators 

such that the form factor each 

b^ depends only on the relative coordinates of particles. 
Any short- range H amiltonian can be explictly expressed 
in this fornP®^®^, which reveals its Laplacian structure. 

Physically, Landau-level-projected Hamiltonians can 
be visualized as a long-range interacting ID chain [Fig.|^. 
The interaction terms can be interpreted as long-range 
(though Gaussian suppressed) hopping processes labeled 
by a. For each a, N particles “hop” from sites rij to sites 
n'-, j = 1,...,A^ according to a COM independent am- 

plitude given hy ^ such that 

the initial and final COM n remains unchanged (Fig.[^). 
Note that although we use the term “hopping”, there is 
no clear distinction between “hopping” and “interaction” 
in our case, as opposed to the Hubbard model. Rather, 
“hopping” is designated for any interaction term that is 
purely quantum (i.e., not of Hartree form). The goal of 
the remainder of this paper is to show how to systemat¬ 
ically construct the polynomial amplitudes Pa that need 
to be inserted into Eq. (10) to obtain the parent Hamil¬ 
tonian for the desired quantum Hall state. 


B. Geometric derivation 

We now specialize Eq. to Wbody interactions that 
are Haldane PPs U^. One appealing feature of the 
second-quantized form of Eq. is that we can construct 
the desired PPs from symmetry principles alone, with¬ 
out referring to the first-quantized form of the interaction 
y(ri,...,rjv). 

The many-body PPs are, by definition, supposed to 
project onto orthogonal subspaces labeled by m, where 
m denotes the relative angular momentum: 

UX = kZ') {^Z'\ IV-^) (V-^l = UZSm'm. (11) 

This requires that 

(0| 6^16™ |0) cx Pm'{Kni,...,KnN) 

+ - • •+'iiiv =0 

/ N \ 

X p„(Kni,...,KnAr)exp 1 

— bm'm- 

If is to vanish with mth total power as N particles 
approach each other, the polynomial Pm must be of de¬ 
gree m. Hence the PPs will be completely determined 
once we find a set of polynomials {pm} such that: (1) 
Pm is of total degree m; (2) Pm has the correct symme¬ 
try property under exchange of particles, i.e. is totally 
(anti)symmetric for bosonic (fermionic) particles; (3) the 


p^’s are orthonormal under the inner product measure 
{Pm',Pm)= E Pm'(Kni, ...,Knjv) 

ni + ...+nAr=0 

Using the barycentric coordinates (Appendix [A[) to rep¬ 
resent the tuple (ni, 

{Pm'.Pm}^ [ dWdQ Pm'{W,Q)pm{W,Q) 

= Sm'm-) 

where W and U are the radial and angular coordinates of 
the vector x G representing the tuple (ni, ...,nAr), 

and Jn -1 is the Jacobian for the spherical coordinates 
in We have exploited the magnetic translation 

symmetry of the problem in quotienting out the COM 
coordinate n. (It is desirable to quotient out n = rij, 
since n takes values on an infinite set when the particles 
lie on the 2D infinite plane, and that complicates the 
definition of the inner product measure.) Each quotient 
space is most elegantly represented as an — 1-simplex 
in barycentric coordinates, where particle permutation 
symmetry (or subgroups of it) is manifest. Explicitly, 
the set of flj can be encoded in the vector 

N 

x = K'^nj^j, (15) 

i=i 

where the N basis vectors {/dj} form a set that spans 
A configuration (ni,...,nAr) is uniquely repre¬ 
sented by a point x G that is independent of 

n = rij. Since x should not favor any particular 

flj, any pair of vectors in the basis {/3j} must form the 
same angle with each other. Specifically, 

~ 

so that each vector points at the angle of tt — cos“^ 
from another. 

With this parametrization, the Gaussian factor reduces 
to the simple form 

AT 

= = (17) 

i=i 

Further mathematical details can be found in the exam¬ 
ples that follow, as well as in Appendix [A| 

The integral approximation in Eq. [^becomes exact in 
the infinite plane limit, and is still very accurate for val¬ 
ues of m where the characteristic inter-particle separation 


X Pm{K,ni,...,K,nN)exp 
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is smaller than the smallest of the two linear dimensions 
of the QH system. In the following, we will assume this 
to be the case; otherwise, there can be significant effects 
from the interaction of a particle with its periodic im¬ 
ages. This was systematically studied in the appendix of 
Ref. [44] for N = 2. We note that the approximation in 
Eq. does not affect the exact zero mode property of 
the trial Hamiltonians constructed below. 


C. Orthogonalization 

We are now ready to evaluate the pseudopotentials. 
To find a second-quantized PPs with relative angular mo¬ 
mentum m, one needs to follow the rules listed in Table [I[ 
In the following, we shall execute this recipe explicitly for 


Write down the allowed “primitive” polynomi- 
ai p0 | 42 | 4^ of degree m <m consistent with the symmetry 
of the particles. 


(2) Orthogonalize this set of primitive polynomials ac¬ 
cording to the inner product measure in Eq. 14 


TABLE I. Summary of the PP construction for spinless par¬ 
ticles. Examples of this procedure are given in Sec. ^ 

= 2, 3, and, to some extent, N = 4-body interactions. 


1. 2-body case 

Eor two-body interactions, we have 

/ oo 

P^^{W)p^{W)e-^^"dW, ( 18 ) 

-OO 

where W = —2nni = 2nn2. Eor this N = 2 case, 
we have allowed W to take negative values as the an¬ 
gular direction spans the ID circle, which consists of 
just two points. The primitive polynomials for bosons 
are { 1 , LE^, LE^, ...} while those for fermions are 
{LE, LE^, LE^,...}. After performing the Gram-Schmidt 
orthogonalization procedure, the 2-body PPs are found 
to be (X where 

^ ^ Pm (^'^1 5 ^ ~ ^nj2-\-ni^nj2-\-n2'> 

ni+^2=0 

( 19 ) 

nl2 ^ rij are integers, and Pm is a mth degree Hermite 
polynomial given in Table |TTj In particular, we recover 
the Laughlin z/ = 1/2 bosonic or 12 = 1/3 fermionic state 
for m = 0 or m = 1 respectively. 

One can easily check that PPs become more delocalized 
in LE-space as m increases. Indeed, 

-I -I noo 

= 0 / Pm{W)W^e---^"dW 

3 ^ ^ 

= TO+1, (20) 


m 

Bosonic Pm (IE) 

Fermionic Pm{W) 

0 

1 

0 

1 

0 

VE 

2 


0 

3 

0 

^{-3 + W^)W 

4 

^{3-6W^ + W^) 

0 

5 

0 

^iu-iow'^+ W'^)W 

6 

^(-15 + 45H/2 

0 


+ IK®) 


7 

0 

^(-105 + lOSlK^ - 211K^ 



+W^)W 


TABLE 11. Representative polynomials pm for the first few 
N = 2-body PPs for bosons and fermions. Pm{W)e~^ 
represents the amplitude that two particles Wjn sites apart 
are involved in a two-body hopping on the chain. 


which is reminiscent of the interpretation of m as the 
angular momentum of a pair of particles in rotationally- 
invariant geometries (Eig. [^. Obtaining the pseudopo¬ 
tentials in this second-quantized form is highly advanta¬ 
geous. In particular, note that this construction is free 
from ambiguities in the choice of E(ri,..., ttv), since sev¬ 
eral possible real-space interactions, e.g., those of the 
Trugman-Kivelson typ^^, can all be grouped into the 
same m sector. This is discussed in more detail in Ap¬ 
pendix 


2. 3-body case 


Eor N = 3, the inner product measure takes the form 

poo p 27 T 

{Pm',Pm)= / Pm'{w,e)p^{w,e)e-i^"wdedw 

Jo Jo 

( 21 ) 

with 

2VE 

Til = cos^ ( 22 a) 

OKj 

VE ^ 

77-2 = -—{Vs sin0 — cosO) ( 22 b) 

3k, 

W r- 

ris = —(—v3 sin^ — cos ^) (22c) 

OK 

Each of the n^’s are treated on equal footing, as one can 
easily check graphically. The above expressions are the 
simplest nontrivial cases of the general expressions for 
barycentric coordinates found in the Appendix (Eqs. |A3[ 


A5). 





















The bosonic primitive polynomials are made up of el¬ 
ementary symmetric polynomials 81 , 82,83 in the vari¬ 
ables hiflj = K{nj — n/N). Since 81 = = 0, the 

only two symmetric primitive polynomials are 

-S 2 = -k }^ninj = - - - = — = -VF , 

i<j i 

(23) 

and 

Y = S3 = K^'[lni 

i 

^3 

= -n2- n3)(2n2 -113- rii) 

(2n3 - ni - 712 ) 

= ^W^cos30. (24) 


The fermionic primitive polynomi als ar e totally anti¬ 
symmetric, and can always be writteii^^^as a symmetric 
polynomial multiplied by the Vandermonde determinant 



FIG. 4. (Color online) The primitive polynomials have beau¬ 
tiful geometric shapes when plotted in V — 1-dim spherical 
coordinates. Shown above are the constant W plots of the an¬ 
tisymmetric polynomial A = (ni — n2)(ni — n3)(ni —n4){n2 — 
^3)(^2 — n 4 ){n 3 — n^) (left), and the degree m = 4 symmetric 
expression in the orthonormalized space spanned by the ele¬ 
mentary symmetric polynomials 1 , 82 , S 3 , 84 and Si (right). 
On the left, there are 24 lobes that each maximally avoid the 
vertices of the tetrahedron (3-simplex). On the right, the 6 
lobes lie around the centers of the 6 edges of the tetrahedron, 
where two of the nds are equal. 


A = K^{ni - n 2 ){n 2 - n 3 ){n 3 - ni) 

= K^(ni - n2)(n2 - n3)(n3 - ni) 

=- ^W^sinSe. (25) 

3v 3 


Note that is of degree 2 while A and Y are of degree 
3. All of them are independent of the COM coordinate 
n, as they should be. N = 3 PPs were derived in Ref. |44| 
through explicit integration, and the approach discussed 
here considerably simplifies those computations by ex¬ 
ploiting symmetry. 

To generate the fermionic (bosonic) PPs up to we 
need to orthogonalize the basis consisting of all possi¬ 
ble (anti)symmetric primitive polynomials up to degree 
m. For instance, the first seven (up to m = 9) 3-body 
fermionic PPs are generated from the primitive basis 
{A, AW^, AY, AW^, AYW^, AY^, AW^}. Note that the 
last two basis elements both contribute to the m = 9 PP 
sector. 

The 3-body PPs are found to be oc 
where 

C= V Pm{W,Y,A)e-"^'^"cn,Cn,Cn„ (26) 

ni+n2+n3=n 


with the polynomials listed in Table m These re¬ 
sults are fully compatible with those from Ref. |42l As 
mentioned, there can be more than one (anti)symmetric 
polynomial of the same degree for sufficiently large m. 
This leads to the degenerate PP subspace, a specific ex¬ 
ample of which is presented in Sec. |V A 2| 


3. N > Yhody case 


For general PPs involving N bodies, the inner product 
measure takes the form 

POO 

{Pm',Pm) = / e-^^"W^-^dW 

Jo 

p2tt N—3 

X / d(pN-2 R / {(j)k)d(j)k 

Jo Jo 

Xpm'{W, 01, ..., (t)N-2)Pm{W, 0i, ..., (j)N-2), 

(27) 


where the Jacobian determinant from Eq. has al¬ 
ready been explicitly included. One transforms the tuple 
{fii, ...,n]sf} into {N — l)-dim spherical coordinates via 
the barycentric coordinates detailed in Appendix [A] 

The bosonic primitive basis is spanned by the el¬ 
ementary symmetric polynomials { 82 , S 3 ,..., Sn} = 


{Hij ^ijk •••, Ylj rij} and combinations 

thereof. For instance, with N = 5 particles at degree 


m = 6, there are 3 possible primitive polynomials: Si, S 3 
and S 2 S 4 . The fermionic primitive basis is spanned 
by all the symmetric polynomials as above, times the 
degree (^) Vandermonde determinant Yli^j{ni — fij) 
shown in Fig. 

From the examples above, one easily deduces the de¬ 
generacy of the PPs to be P{m,N) — P{m — 1 ,N) 

for bosons and P (^m — (^), — P (^m — 1 — (^), vj 


for fermions, where P{m,N) is the nu mber of partitions 
of the integer N into at most m partJ^^EH. In particu¬ 
lar, the degeneracy is always nontrivial (> 2) whenever 
V > 4 and m > 4. 
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m 

Bosonic pm{W,Y) 

Fermionic pm{W, Y, A) 

0 

1 

0 

1 

0 

0 

2 

i-lw^ 

0 

3 

3V2Y 

A 

4 

1 - 

0 

5 

V2Y{6 - W^) 

A (2- IW^) 

6 

(i) 1 - 2W'^ + 

(ii) 7291^2) 


7 

^{45 - 15W^ + W*)Y 

VlOA (1 - ^^4^ 

8 

(i) 1 + - 6)(54 - 18W'^ + 

(ii) 2i3)/¥(2^^ “ 21)(2VF® - 729^^) 

1 

CO 

9 

(i) y^y(-30 + i5W^ - 21^^ + isr^) 

(ii) ^y^y(-5760VF2 + 7567^^ - 31W^ + 81(140 + 9^^)) 

(i) ^i-10 + 5W^-IW^ + iW^) 

729y^) 


TABLE III. The polynomials pm for N = 3-body PPs for bosons and fermions up to m = 9. Note that there is more than one 
possible PP for larger m, since there are multiple ways to build a homogeneous (anti)symmetric polynomial from the primitive 
and elementary symmetric polynomials. For instance, with Y and A dehned in Eqs. |24| and there are two possible bosonic 
PPs for m = 6, since there are two ways {W^ and to build a homogeneous 6-degree polynomial from elementary symmetric 
polynomials. For the fermionic cases, we have explicitly kept only one factor of A, since even powers of A can be expressed in 
terms of and Y. For all cases, the mean-square spread of the PPs also increases linearly with m, i.e. ^‘j) = + 1 

(compare with the two-body case in Eq. 20). 


IV. GEOMETRIC CONSTRUCTION OF 
PSEUDOPOTENTIALS: SPINFUL CASE 

In the presence of “internal degrees of freedom” 
(DOFs) which we also refer to as “spins” or “compo¬ 
nents” for simplicity, there are considerably more pos¬ 
sibilities for the diverse forms of the PPs. This is be¬ 
cause the PPs consist of products of spatial and spin 
parts, and either part can have many possible symme¬ 
try types, as long as they conspire to produce an overall 
(anti)symmetric PP in the case of (fermions) bosons. 

A generic multicomponent PP takes the form 

( 28 ) 

A cr 

The notation here requires some explanation. A = 
[Ai,A 2 ,...] defines a partition of all particles into sev¬ 
eral subgroups Xi each of which imposing a symmetry 
constraint among particles of this subgroup. Associated 
with it is which is a spatial term (in (ni,..., nTv)- 
space) exhibiting this symmetry. |cr^) = |cricr 2 ...) refers 
to an (internal) spin basis that is consistent with the 
symmetry type A. Each symmetry type corresponds to a 


partition of V, with Xj = N. For N bosons, A repre¬ 
sents the situation where there is permutation symmetry 
among the first Ai particles, among the next A 2 parti¬ 
cles, etc. but no additional symmetry between the Xi 
subsets. This is often represented by the Young Tableau 
with Xi boxes in the i^^ row. For fermions, we use the 
conjugate representation A, with the rows replaced by 
columns and symmetry conditions replaced by antisym¬ 
metry ones. For instance, the totally (anti)symmetric 
types are [1,1,..., 1] and [Y], respectively. 

Hence the space of PPs is specified by three parame¬ 
ters: N - the number of particles interacting with each 
other, m - the total relative angular momentum or the 
total polynomial degree in and b - the number of inter¬ 

nal DOFs (spin). While the symmetry type A and hence 
^ depends only on N and m, the set of possible |cr^) 
also depends on b. To further illustrate our notation, 
we specify N^m^b parameters for the interactions rele¬ 
vant to some commonly known states: in the archetypical 
single-layer FQH states, we have 6 = 1 components, and 
the N = 2 PP interactions for the Laughlin state penal¬ 
ize pairs of particles with relative angular momentum m, 
where is the filling fraction. For bilayer FQH states, 
we have 6 = 2 and N = 2-body interactions. PPs as en- 
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ergy penalties in the [2] sector with angular momentum 
< m and [1,1] sectors with angular momentum < n give 
rise to the Halperin {mmn) states. Here, the [2] sector 
is also known as the triplet channel, as it is spanned by 
the following three basis vectors: {|tt) : Iti) + lit) : lit)}- 
By contrast, the [1,1] sector only contains {|tt) ~ lit)} 
as dictated by antisymmetry. 


A. Multicomponent pseudopotentials for a given 
symmetry 

The construction of multicomponent PPs here paral¬ 
lels that of multicomponent wave functions described in 
Ref. [43l For completeness, we first review this construc¬ 
tion, and proceed to show how an orthonormal multi- 
component PP basis, adapted to the cylinder or torus, 
can be explicitly found through the geometric approach. 
We describe how to first find the spatial part of the PP 
and second the spin basis |cr^)- 


1. Spatial part 

For each symmetry type A, we can construct the spa¬ 
tial part ^ with elementary symmetric polynomials in 
subsets of the particle indices fii = rii — n/N. They are, 
for instance, 5'ip2 = hi+h2, 5'2,234 = h2h3-hh3h4+h2h4, 
etc. Of course, we must have S'!,123. ..at = = 0. 

Like in the single-component case, the spatial part 
consists of a primitive polynomial which enforces 
the symmetry, and a totally symmetric factor that does 
not change the symmetry. Here, the main step in 
the multicomponent generalization is the replacement of 
primitive polynomials 1 and A by primitive polynomi¬ 
als consistent with the symmetry type A. As the sim¬ 
plest example, the primitive polynomial in [2 1] 

Si ^12 = —fis (and cyclic permutations). It is the only 
possible degree m = 1 expression symmetric in two (but 
not all three) of the indices. 

In general, there can be more than one candidate 
monomial obeying a symmetry consistent with A. For 
instance, for U^^ 2[2 1] ^1,12 = (^1 + ^2)^ and 

5^2,12 = nin 2 (and cyclic permutations thereof). To find 
the primitive polynomials, we will have to construct one 
or more linear combinations of these terms which do not 
have any higher symmetry other than [2,1] (i.e. in this 
case, this higher symmetry channel could be [3]). El¬ 
ementary computation reveals that the only primitive 
polynomial should be Sf 12 + ‘^^ 2 , 12 ^ because it is the 
only linear combination that is manifestly symmetric in 
indices 1, 2 and disappears upon symmetrization over all 
three particles. 

A = 4, m = 2, A = [2, 2]. As a more involved exam¬ 
ple, we demonstrate how to find the primitive polynomial 
corresponding to the symmetry type [2,2]. Independent 
monomials that satisfy this symmetry include Sf 12 ^ 5'2,i2 


and 5'2,34- The primitive polynomial is then given by the 
linear combination 


^1,12 + l ^ lS 2,12 + /^25'2,34 (29) 

with fii and fi 2 to be determined by demanding that the 
linear combination disappears upon symmetrizing over 
permutations under [3,1] and [4]. The symmetrized sums 
are 

2((^2 + '^3)^ T (^2 + '^4)^ + (^3 T '^4)^ T 
/ii(nin2 + nins + 712^3) + /a 2 {nin 4 -h 712^4 + 72377.4)) 

(30) 

and 


4 (722 + 723 + 72I + 722723 + 723724 -h 722724) (4 - /2i - 1 ^ 2 ). 

(31) 

Setting them both to zero, we find /ii = /i2 = 2, so the 
primitive polynomial is Sf 12 + 25'2,i2 + 25'2,34- 

A = 3; 772 = 2, A = [2,1]. The above procedure 
works for arbitrarily complicated cases, but quickly be¬ 
comes cumbersome. This is when our geometric approach 
again becomes useful. We first write down the relevant 
monomials in barycentric coordinates given by Eq.j^for 
A = 3 bodies (or Eq. A12 for general A). The coefficients 
in the primitive polynomial can then be elegantly de¬ 
termined through graphical inspection. We demonstrate 
this explicitly by revisiting the example on A = [2,1]. 
Recall that the primitive polynomial (call it Prn =2 = P 2 ) 
is a linear combination of Sf 12 ^2,12, i-e. 

P2 = ^1,12 T /^^ 2,12 


= — 2 + (1 + /2)(cos 26> — A/3sin26>)^ 


(32) 


where 0 = 0,27r/3,47r/3 points towards the vertices fa¬ 
voring 72i, 722,723 respectively. The correct value of /i will 
cause P 2 to disappear under symmetrization of the 3 par¬ 
ticles. Graphically, it means that the lobes of three copies 
of the plot of ^ 2 ^ each rotated an angle 27r/3 from each 
other, must cancel upon addition. This is illustrated in 
Fig. H where /i = 2 is readily identified as the correct 
value. 

All in all, we have the 772 = 2 primitive polynomial for 

[ 2 , 1 ]: 


P 2 = (f^3)^ + 272ih2 = COS 20 + sin 20j . (33) 


A = 3; 772 = 1, A = [2,1]. The reader is invited to also 
visualize the simpler case for the [2,1] sector, which takes 
the form 


Pi = Si^i2 = —723 


w 

3k 


^cos 0 -f 



(34) 


Compared to the antisymmetric case for fermions, we 
have two primitive polynomials P 2 , instead of just A. 




11 



FIG. 5. (Color online) Polar plots of P2{0) as a function of 0^ 
with /i = 0.5 (blue curve) and /i = 2 (purple curve). ^2(0) + 
^2{0 + 27r/3) + yd 2 (^ + 47r/3) sum to zero only when the IoIdcs 
are of equal size, which is the case for /x = 2 only. 

The primitive polynomials for a wide variety of cases are 
also listed in Table II and III of Ref. [43l and can alterna¬ 
tively be systematically derived via the theory of Matric 
units explained in the appendix of the same reference. 

2. Admissible spin bases 

In general, not all possible A^-particle spin bases sur¬ 
vive under symmetrization with respect to a given sym¬ 
metry type A. Only those that survive should be in¬ 
cluded, i.e. are admissible, in the set of basis in 
|cr^). For instance, the basis \aaa) does not appear 
in To see this, try symmetrizing the combi¬ 

nation p(ni,n 2 ,%) \aaa) subject to the condition that 
p(ni,n 2 ,^ 3 ) has no higher symmetry than A = [ 2 , 1 ], i.e. 
the totally symmetrized sum S [p(f^i,^ 3 )] = 0- Obvi¬ 
ously, p(ni, n 2 , '^ 3 ) \aaa) then has to symmetrize to zero 
and should not be included in [ 2 , 1 ]. 

There is a nice way to write down the set of admissible 
bases by looking at the labelings of semistandard Young 
Tableaux. From Ref. sa the admissible bases in |cr^), up 
to permutations, are in one-to-one correspondence with 
the labelings of semistandard Young Tableaux with num¬ 
bers 1 to N. (A semistandard Young Tableau has non¬ 
decreasing entries along each row and strictly increasing 
entries along each column.) For instance, a symmetry 
type of [3,1] with 6 = 3 corresponds to labelings (with 
labels a, 6 , c defined to be in increasing order): 

[aaa, 6 ], [aa 6 , 6 ], [aac, 6 ], [a 66 , 6 ], [a 6 c, 6 ], 

[acc, 6 ], [aaa, c] , [aa 6 , c] , [aac, c] , [a 66 , c ], 

[a 6 c, c] , [acc, c], [ 666 , c], [ 66 c, c], [ 6 cc, c] and [ccc, c], 

where rows are separated by commas. From these, the 


admissible states for A = [3,1] can be written down by 
copying the labelings verbatim: 

|aaa6), |aa66), |aac6), |a666), |a6c6), 

|acc6), |aaac), |aa6c), |aacc), |a66c), 

|a6cc), |accc), |666c), |66cc), |6ccc) and |cccc). 

In the Young Tableau corresponding to A, each row rep¬ 
resents a set of particles that are symmetric with each 
other in the A representation. The requirement that la¬ 
bels increase monotonically within each row defines an 
ordering, and prevents the repeated listing of basis states 
related by permutation. The requirement of strictly in¬ 
creasing labelings down each column also prevents that, 
and also avoids the listing of bases that do not survive 
under symmetry constraints. 

With these preliminary considerations, we are now in 
position to formulate the general recipe for constructing 
PPs in the multicomponent case, which is given in Ta¬ 
ble |IV| In the following section we apply this recipe to 
the case of SU(2) spins with 3-body interactions. 


(1) Given N,m and 6, specify which symmetry type the 
PP is associated with. For example, we specify the [2] 
or the [1,1] channel when computing PPs realizing the 
fermionic Halperin states in bilayer QH systems (6 = 
2, A = 2). 


(2) Next, determine the appropriate primitive polyno¬ 
mials by hnding the coefficients multiplying the allowed 
monomials. Then multiply the primitive polynomial by 
symmetric polynomials, and orthogonalize to obtain the 
spatial parts U^ x- 

(3) Finally, choose the desired admissible spin channels, 
and (anti)symmetrize the resultant product of the spatial 
and spin parts depending on whether we want a bosonic 
(fermionic) PP. 


TABLE IV. Summary of the PP construction procedure for 
the multicomponent case. Examples of this procedure are 
given in Section |V| 


B. Example: N = 3-body case for SU(2) spins 

Here we explicitly work out the multi-component case 
6 = 2 which is also the most common multi-component 
scenario. We focus on A = 3-body interactions to illus¬ 
trate the nontrivial aspects of our approach. 

The symmetry type is given by [Ai,A 2 ] = [Y/2 + 
S', A/2 — S], where S is the total spin of the particles. 
Let us denote the spins by There are only A = 3 
boxes in the Young Tableaux, with the following possible 
symmetry types: [3], [2,1], [1,1,1]. 

For type A = [3], the possible spin bases are Ittt)? 

Itii) and For A = [2,1], the possible 
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spin bases are Itii), Itti)? corresponding to Tableau la¬ 
bellings and respectively. For A = [1,1,1], 

there is actually no admissible spin basis: total (internal 
DOF) antisymmetry is impossible for 3 particles, if there 
are only 2 different spin states to choose from. 

Consider bosonic particles in the following. The sym¬ 
metry type [3] case corresponds to the primitive polyno¬ 
mial 1, so the resultant PP takes the same form as in 
the single-component case. After symmetrizing the spin 
part, we have the following available sets of spin channels 
for [3]: 

{ittt)}, 

{im), int), int)}, 

For the PP pj oc we have 

c^|0)= E (35) 

for S = 3/2, and 


m 

Pm,[2,l]{Pl,p2,W,Y) 

1 

I3i 

2 


3 

{-3 + W^) 

4(i) 

^(/?2 + 6/?ir) 

4(ii) 

^ (-54^2 + 12 I32W^ + (cos[46»] -h ^/3sin[46»])) 

5(i) 

^(/3i (-3 + 2W^)+6/32Y) 

5(ii) 

cos[56'] -h 3(5^/3/3l(27 - 181F^ 


-h2VF^)-1-IF® sin[56»])) 


TABLE V. The polynomials pm for the first few N = 3-body 
PPs for bosons in the total spin \S\ = |, i.e. A = [2,1] chan¬ 
nel. The primitive polynomials /3i and p 2 are given by Eqs. |33| 
and [34] cos 4^, sin 4^, cos 5^, sin 5^ can all be decomposed into 
the elementary symmetric polynomials and primitive polyno¬ 
mials IT, y, /3i, and ^2 ■ 


b]r\o) = s 


p^e im) 


(36) 


for S' = 1/2. Other contributions with S ^ —S can be 
obtained via the identification |t) ^ |i). Here |tti) is 
the shorthand for |0)> P™ = Pm,[3] 

refers to the same polynomial as in the single-component 
case. 

To construct the orthonormal PP basis for 
the first few m, we orthogonalize the set 
{/?l,/32,/?lVF^/32Tr^/3lr,/32r,/3lTP^...}. The re- 
suits are shown in Table |Vl Notice that the first PP 
for symmetry type [2,1] occurs at m = 1, whereas that 
of single-component bosons/fermions occur at m = 0 
and m = 3 respectively. Indeed, there are more ways of 
constructing PP polynomials when only a subgroup of 
the full symmetric/alternating group is involved. The 
onset of PPs degenerate in m also occurs earlier, at 
m = 4. 

To illustrate the full procedure for PP construction 
with internal DOFs, we detail the case of S = —1/2 
below. For m = 1, 




Y (-h 3 )e-^'^^m) 


= -e-3'^'[(h2 + h3)|m) 

-h {fii + ns) liti) + + ^2) liit)] 

= e 3 ^ [ni Itii) + n 2 |iti) + % lilt)]? 


( 37 ) 


while for m = 2 we have 


6^t|o)=5 


^ ((h3)2 + 2hih2)e-5'^'|m) 


T.jnj=n 

= e“^^^[((n2)^ + (ns)^ + 2ni(n2 +%)) |tU) 

+ ((f^i)^ + {^ 2 )^ + 2n3(ni + n2)) |iit) 

+ ((f^i)^ + (i^s)^ + 2n2{ni -h ns)) liti)]: 


(38) 


Expressions for S' = 1/2 are obtained via |t) ^ |i). 

For the present case of SU(2) spins, i.e. 6 = 2, there 
exists a nice closed-form generating function for the di¬ 
mension of the spatial basis for each m. Define the gen¬ 
erating function = Y.m d{[\i,\ 2 ]]m)q^ where 

c^([Ai, A 2 ]; m) is the number of different polynomials with 
degree m and symmetry type [Ai,A 2 ]. It can be shown 
that for bosons. 


^[Ai,A2](9) 


__ 

nEi(i-9'")nE(i-3") 
_ 1 - g _ 

nEi(i-9’”)nAAi-</”) 


(39) 


This is obtainecP^by considering the dimensionality from 
two symmetry subsets, and then subtracting overlaps 
from the higher symmetry case [Ai + 1,A2 — 1]. The 
dimension is related to the q-binomial coefficient. For 
6 > 2, however, the situation is much more complicated, 
involving Kostka coefficients which count the number of 
semisimple labelings of A with a given alphabet. 
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V. PSEUDOPOTENTIAL HAMILTONIANS: 
CASE STUDIES 


generalized Moore-Read Pfafiian state at filling v = \jq^ 
for which the g = 4 case will be discussed in detail. 


In this Section we provide several applications of the 
general pseudopotential construction developed in the 
previous sections, with examples arranged in the order 
of increasing com plexity. We start with spin-polarized 
states (Sec. VA), whose Hamiltonians have been ob¬ 
tained via alternative methods and are well-known in the 
literature. As an illustration of the method, we provide 
a detailed derivation of the fermionic Gaffnian parent 
Hamiltonian. Note that the resulting second-quantized 
form of the Hamiltonian applies (with minimal modifi¬ 
cations) to both cylinder and torus geometries. In ad¬ 
dition to pedagogical examples which illustrate our ap¬ 
proach, we construct pseudopotential Hamiltonians for 
some non-Abelian states that have not been available in 
the literature for any geometry. 

The main sequence of steps is stated as follows: 


1. Assuming that the ground-state wave function is 
known in the first-quantized form, find its thin 
torus pattern (this step has been frequently dis¬ 
cussed in the literature; for completeness, we pro¬ 
vide a brief summary in Appendix [C|). 

2. Using the formalism developed in the previous 
sections, write down the parent Hamiltonian cor¬ 
responding to the root pattern in the second- 
quantized form. 


3. Use numerics (exact diagonalization) to verify that 
the ground state of the proposed Hamiltonian is 
indeed given by the initial wave function. The ver¬ 
ification criteria include testing for the unique zero- 
energy ground state at the given filling factor, the 
correct thin-cylinder root patterns as the circum¬ 
ference of the cylinder is taken to zero, and the 
level counting of entanglement spectra. For the 
purpose of numerical calculations, we place a fi¬ 
nite total number of particles N on the surface of 
a torus or an open cylinder. The two linear di¬ 
mensions of the Hall surface {L and H) satisfy the 
relation LH = 27r/^Aorb: where Aorb is the number 
of available orbitals (it is equal to on the torus, 
and equal to A^ -h I on the cylinder). Unless stated 
otherwise, we also assume L = H. 


A. Spin-polarized states 

The construction of two-body as well as the shortest- 
range three-body pseudopotential Hamiltonians has been 
discussed in depth in the literature, see e.g.. Refs. 01 El- 
[TTl Starting from there, as our most elementary example 
for spin-polarized particles, we consider ground states of 
longer ranged 3-body potentials: the fermionic Gaffnian 
state at filling u = 2/5 as well as the Haffnian and the 


1 . Fermionic Gaffnian 


The derivation of the fermionic Gaffnian parent Hamil¬ 
tonian is summarized in Fig.[^ The wave-function of the 
Gaffnian state on the infinite plane is given 

Vl/Gaf = .4{^'332({^t}, {Zi})Per }, (40) 

where we have suppressed the usual Gaussian factors. 
Although physically a one-component (spin-polarized) 
state, the above wave-function reflects the underlying 
“two-component” nature of the Gaffnian state: in order 
to write down the wave function, we have divided elec¬ 
trons into two groups t and | (not to be confused with 
physical spin) that are correlated through the Jastrow 
factors within the 332 Halperin state, which is defined as 


«'332({2:t> 2^4.}) = n(^t,i - Zfjf n(^-^-» “ 

i<j i<j 

(41) 

hj 

as well as the “permanent”, 


Per 


= y _^_. 

zt-ziJ yYlMj-ziau)) 


(42) 


Ultimately, the distinction between ^ and ^ particles is 
erased by the overall antisymmetrization M, producing a 
well-defined single-component wave function. 

Having obtained the first quantized wave function for 
the Gaffnian, the second step is to determine its thin 
cylinder root patterns. (Note that a similar, but dif¬ 
ferent analysis can also be executed on the sphere by 
analyzing the root partitions of the parent Hamiltonian 
null spac^^ . The detailed procedure for finding the root 
patterns of Halperin bilayer states was given in Ref. [56l 
of which a brief summary is outlined in Appendix The 
Gaffnian wave function vanishes as power 6 as three par¬ 
ticles are brought together, hence its thin torus root pat¬ 
terns are IIOOOIIOOO... and lOIOOIOIOO..., accordingly. 

From the form of the root patterns, we conclude that 
we need and terms in order to build the par¬ 

ent Hamiltonian. The polynomial amplitudes for these 
terms can be found in Eq. [^and Table [Hlj The Gaffnian 
parent Hamiltonian therefore reads 


H = ybibn+-fydidn, ( 43 ) 

n n 

E /^^ — 2 

A{ni,n2, nz)e~^^i=^ (44) 

ni+n2+n3=0 


dn= E A{ni,n2,ns) 

ni+n2+n3=0 

-2 

xe 2 = 1 JCn^Cn^Cns, 


2 - gW^^(ni,n2,n3) 


(45) 
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where as before fij = rij — n/3. The Hamiltonian (43) 
assigns positive energies to m = 3, 5 in a cluster oF3 
particles in an infinite system, while all the other ener¬ 
gies are exactly zero. The constant 7 tunes the ratio of 
the two non-zero energies and can be set to any posi¬ 
tive number. The precise value of 7 does not affect the 
ground state wave function or its energy, but it does have 
an effect on the energetics of the low-lying excited states. 
In the following, for the sake of brevity, we refer to the 
Hamiltonian in Eq. (43) simply by H + 7^(2 — 


Finally, having obtained the parent Hamiltonian, we 
can verify that it yields the correct ground state we 
started from (Eq. 40). As we see in Fig. exact di- 
agonalization on the torus for = 6 , 8,10 electrons con¬ 
sistently finds a zero-energy ground state that is 10 -fold 
degenerate and has zero momentum (i.e., does not break 
translation symmetry). One can further verify that we 
have obtained the correct ground state by computing its 
entanglement properties, as we discuss below. Another 
non-trivial check is to perform diagonalization on a torus 
stretched along one axis, and directly identify the root 
patterns in the thin-torus limit. In doing so we find, 
as expected, the ground states to evolve to the single 
Fock states 1100011000... and 1010010100... (and those 
obtained from these by COM translation), in agreement 
with the fermionic Gaffnian state. 


(1) First quantized wave function: 

(4) Energy spectrunn on the torus 

^Gaf = ^{^332per^^:4^} 

0 . 7 1 


1 

0.6 ■ 

(2) Thin-torus root patterns: 

...110001100011000... 

...101001010010100... 

0.5 ■ 

LU 0.4 

0.3 



0.2 

- 

(3) Parent Hanniltonian: 

0.1 

10-fold degenerate - 

A,A{2-^W^) 

0 ' 

ground states 

( 

D 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1 


klg 


FIG. 6 . (Color online) Main steps in constructing the 
fermionic Gaffnian parent Hamiltonian, (i) Assume the first 
quantized wave function, in this example given by Eq. |40[ 
(ii) Find the corresponding thin torus root patterns (here 
1100011000... and 1010010100...). (hi) Construct the par¬ 
ent Hamiltonian by consulting Eq. (26) and Table |HI| (iv) 
Verify the construction and compute properties by exact di¬ 
agonalization. Here, this is performed on a hexagonal torus 
with N = 6,8,10 particles, finding a 10-fold ground state 
with exactly zero energy. The energy scale is set by 7 = 1 
in Eq. As the torus is stretched along one direction, the 
ground states evolve to the expected patterns 1100011000 ... 
and 1010010100.... In addition, we obtain the entire low-lying 
neutral excitation spectrum, allowing one, in principle, to ex¬ 
tract many other properties of the system. 


2. Pfaffians and Haffnians 


We next conside r the generalized non-Abelian Pfafhan 
state at 1/q fiihn^Hini. The wave function in the disk 
geometry reads 

where even (odd) q corresponds to a fermionic (bosonic) 
state. The Pfafhan is dehned as 

^ n/2 

^ 2»/2(n/2)! ^ sgn(a)n^<x( 2 i-i),a( 2 i), 


where A is the nxn skew-symmetric matrix Aij = 1 /( 2 ;^ — 
Zj) for n even. The case q = 2 reduces to the familiar 
Moore-Read state, which is the ground state of purely 
3-body interaction l/f. One would naively expect that 
q> 2 states can be obtained by adding 2 -body terms to 
t/ 3 , but as one can explicitly verify, this is not the case. 

By a power counting procedure^ (further elaborated 
on in Appendix the root conhguration of the 1/q- 
Pfafhan state is given by 

10^“^10n0^“^10^. 10 ^“^ 10 T.. (47) 


where 0^ represents a string of q zeros. It follows that its 
parent Hamiltonian is given by the 3-body PP 

^m=3(g-l)’ 

as well as all nonzero spinless 2-body PPs 

K<,-2- (49) 


For the more interesting g = 4 case which we have studied 
numerically, we require the two-body Uf PP and three- 
body and PPs given by Eq. 26 and Table 

|HI| (Again, we emphasize that t/’s promotedto how they 
appear in the Hamiltonian should be understood as poly¬ 
nomial amplitudes that enter the definition of operators 
bl,bn in Eq.[lol) 

Note that labels (i),(ii) stand for two linearly inde¬ 
pendent 3-body PPs that occur for m = 9. Thus, the 
1/4 Pfafhan state is an example whose parent Hamilto¬ 
nian contains a degenerate subspace of PPs. The state 
can alternatively be obtained numerically through Jack 
polynomials^ and the root conhguration given above. 
We have conhrmed that the overlap between the ground 
state of the parent Hamiltonian and the Jack polyno¬ 
mial wave function is equal to one (within machine pre¬ 
cision) for q = 4 and N = 6,8 electrons. Note that 
experiments^ hnd some evidence for an incompressible, 
potentially non-Abelian, = 1/4 state. This may be at¬ 
tributed to the Pfafhan 1/4 state, although theoretical 
calculations suggest that the Halperin 553 state is also a 
candidatS^. 
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For our final spin-polarized example we consider the 
fermionic Haffnian stat^^, whose bosonic counterpart 
was recently considered in Ref. [39l The fermionic 
Haffnian occurs at the filling fraction v = 1/3, and, up 
to COM translation, has the following root patterns^ on 
the torus: 110000110000... and 100100100100.... The lat¬ 
ter pattern is identical to that of the Laughlin state. The 
fermionic Haffnian vanishes as power 7 as three particles 
are brought together, hence we need to impose an en¬ 
ergy penalty for m = 6 too, and the parent Hamiltonian 
consists of 


f/|, Vl, and 


which are given by Eq. 26 and Table HI 


(50) 


B. Spinful states 

We have previously mentioned that if we allow the spin 
degree of freedom to enter, the number of possible states 
obviously becomes much richer with interactions involv¬ 
ing three or more particles. A systematic investigation of 
these states is left for future work. Here we content our¬ 
selves with illustrating our method by formulating parent 
Hamiltonians for several states that have been the sub¬ 
ject of recent attention: the Spinful Gaffniarl^ and the 
non-Abelian Spin Singlet (NASS) stated. The Hamilto¬ 
nians for these states have previously been written down 
for the sphere (or disk) geometry, which we now extend 
to the cylinder and torus. Furthermore, we propose the 
parent Hamiltonian for a certain type of state involv¬ 
ing the permanent (“221 times permanent” stat^^, for 
which the Hamiltonian was previously unknown. 


and the total number of orbitals I a in the left subsystem 
remain good quantum numbers. On the other hand, for 
the PES, the translation or rotation symmetry of the full 
system is preserved. For an open cylinder, as in Fig. 
this means that the total momentum Ky of the left sub¬ 
system is also a good quantum number: 

Ky = '^ mhm, (51) 

m^A 


with fim the density operator acting on the momentum 
orbital m. This linear momentum on the cylinder cor¬ 
responds to the Lz projection of angular momentum on 
the sphere. The PES on the sphere, however, is also in¬ 
variant under full SU(2) rotation, i.e. it forms multiplets 
of L^. This symmetry is absent on a cylinder, where we 
can only use Ky to classify the levels of the entanglement 
spectrum. This is discussed in more detail in the exam¬ 
ple in Sec. VB 1 As we perform an orbital partition, 
the symmetry of the subsystem is reduced, and even on 
the sphere, the only remaining quantum number is Lz. 
Thus, the OES on the sphere can be directly compared 
with the OES on the open cylinder. Finally, for spinful 
states, the spin quantum number commutes with the re¬ 
duced density operator of the OES, and thus allows for 
additional resolution of the ES level counting^. 


1 . Spin-singlet Gaffnian state 

The spin-singlet Gaffnian state is a nontrivial spin¬ 
ful generalization of the bosonic spin-polarized Gaffnian 
state by Davenport et alW^. In addition to the two spin 
polarized 3-body terms with S = 3/2, its parent Hamil- 



Na,Ia,K^ 



a 


FIG. 7. (Color online) Two choices of partitioning the sys¬ 
tem for computing the entanglement spectrum. Orbital cut 
(left) amounts to partitioning the system into two groups of 
orbitals A and B, and tracing out the orbitals in B (marked 
yellow). The Schmidt levels of A can be classified by the num¬ 
ber of particles Na, the number of orbitals I a, and the total 
momentum Ky. An alternative choice is particle partitioning 
(right), where some particles (marked in yellow) are traced 
out, regardless of their position. 


In addition to formulating the parent Hamiltonians, 
we also computed the orbital (OES) and particle entan¬ 
glement spectra (PES) for the respective states. Fig. 
illustrates the two types of partitioning in the case of an 
open cylinder. After performing the cut, the entangle¬ 
ment spectrum can still have some remaining symmetry 
which can be used to classify the Schmidt levels. For 
example, for the OES, the total number of particles Na 
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FIG. 8. (Color online) PES of the spin-singlet Gaffnian state 
for A = 12 particles and = 12 flux quanta on an open 
cylinder with aspect ratio equal to 1. The spectrum is ob¬ 
tained by tracing out Nb — Na — 6 particles and plotted as 
a function of momentum Ky in part A. Different values of 
the total z component of spin in part A, are indicated in 
the inset. 
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tonian also contains the shortest range (m = 1) term 
with S = 1/2. In total, this encompasses the PPs 


^0,5=3/2 5 ^2,5=3/2 ^1,5=1/2- (^2) 


As such, the spin-singlet Gaffnian wave function vanishes 
as the third power in the S' = 3/2 channel, and the sec¬ 
ond power in S = 1 /2. These projectors are given by 


Eqs. 


35 


and Table 


III 


for I/q 3/2 ^2 3 / 2 ’ 


37 


for U/s=i /2 (with both spin orientations |t) ^ |i))- In 
addition to these, we also add the total spin operator 
to our Hamiltonian, to ensure that the ground state is a 
spin-singlet. 


By diagonalizing the Hamiltonian (52) numerically, we 


find a zero-energy ground state at filling factor zy = 4/5 
and shift of —3 on a finite cylinder for N = 4,8, and 
12 particles. We have furthermore computed the PE^ 
for N = 12 particles and = 12 flux quanta, shown 
in Fig. As illustrated in Fig. (right), to perform a 
PES “partition”, we divide the system into parts A and 
H, which both contain orbitals, but Na and Nb par¬ 
ticles, respectively, such that Na + Nb = N. To obtain 
Fig-i we have traced out Nb = 6 particles from the sys¬ 
tem. The PES obtained in this way provides information 
about the counting of quasihole excitations of the given 
state, as shown in Ref. nn 


We compare the counting in Fig. with the corre¬ 
sponding PES obtained on the sphere in Davenport et 
aM Note that our result in Fig. superficially looks 
different from the result in Ref. [62l This is because 
the PES partition preserves the symmetry of the ground 
state, causing the PES on the sphere to have exact ro¬ 
tational symmetry which is absent for our open cylinder. 
In other words, on an open cylinder, the good quantum 
number after partitioning is only the linear momentum 
K^. This quantum number, in turn, corresponds to the 
Lz projection of angular momentum on the sphere. This, 
however, does not exhaust all symmtries of the PES on 
the sphere, where the full angular momentum is a 
good quantum number. This additional degeneracy of 
the PES is factored out in Ref. |62| Furthermore, on both 
the sphere and cylinder, due to the singlet property of 
the wave function, we ex pec t the PES levels to be mul- 
tiplets of the 5'^ operator^, as can be verified in Fig. 
The most important universal information is the count¬ 
ing of PES levels per momentum sector, which we find 
to be in agreement with Ref. [62] As we restore the 
degeneracy in the PES given in Ref. [62] the counting for 
the first three sectors is a single level with S = 1 and 
= 22, 4 levels at = 21 (two with S = 1 and two 
with S' = 0), and 10 levels at Lz = 20 (one of them with 
S = 2, 6 with S = 1 and 3 with S = 0). While the finite- 
size splitting between these levels is non-universal (and 
differs between sphere and cylinder), we indeed obtain 
the identical counting per sector (Fig. [^. 


2. The NASS state 

Another class of non-Abelian spin-singlet states has 
been proposed by Ardonne et under the name “non- 
Abelian spin singlet” states (NASS). The bosonic family 
of such states occurs at filling factors u = 2k/3. Ac¬ 
cording to Ref. |64l the NASS wave function ought to 
vanish quadratically when we bring together k 1 par¬ 
ticles of the same spin, and linearly for /c + 1 particles 
of different spins. For SU(2) spins, such wave func¬ 
tions only exist in the [k + 1] and [/c, 1] representations. 
Its parent Hamiltonian should therefore contain the PPs 
C/*=o,S=fe /2 and f7^=o,S=fe/2-i- since the simplest non- 
trivial totally/partially symmetric polynomial are of de¬ 
grees two and one, respectively. Explicitly, the 
operators are 

^n!s=fe/2 |0) = 6“"^ in - t) , (53) 

^ 5 =^/ 2 -! io) = e-^ ^'”?[in - n) + in - n) 

+ in-nt) + - + in-t)], ( 54 ) 

and similarly for the terms S ^ —S obtained by ex- 
changing |t) |4-)- 

The simplest member k = 2 of the NASS family has re¬ 
cently been studied from the “squeezing” perspective^. 
Below we complement these results by independently 
generating the NASS k = 2 state using its parent Hamil¬ 
tonian in Eqs. [53]and|54| We compute its PES and OES 
on a finite cylinder. Fig. The PES and OES are com¬ 
puted for A^ = 12 bosons on a cylinder with aspect ratio 
1. Similar to the previous cases, the universal informa¬ 
tion in the PES spectrum - the counting of levels per 
momentum sector - is in agreement with the indepen¬ 
dently obtained data in Ref. [64] on the sphere. In this 
case, the bipartition of the system is obtained by fixing 
the cut such that there are /a = 4 orbitals in part A. We 
compute the OES for two choices of the total number of 
particles in A, Na = 6 and Na = 7. We see that the 
counting of the OES changes as a function of Na^ as it 
often happens for non-Abelian states where moving the 
cut probes different topological sectors of the theory, in 
this case associated with SU(3)2 GET. 


3. Halperin-permanent states 


Finally, we tackle more complex examples of non- 
Abelian states that involve a product of Halperin mul¬ 
ticomponent state^^and the permanent. To be specific, 
we study states of the following form in the disk geome¬ 
try: 


4 > 


lln 


Per 




(55) 


where Zi^Wi refer to the spin up and spin down particle 
coordinates, respectively, and i/un is the Halperin wave 
functioiP^, an example of which was given in Eq. for 
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FIG. 9. (Color Online) PES (left) and OES (right) for the NASS z/ = 4/3 state. The system contains A" = 12 bosons on a 
finite cylinder with aspect ratio 1. Eor PES, the subsystem A contains Na = 6 particles, and PES is plotted as a function of 
momentum Ky of subsystem A. The states are labeled by different values of S^, the total z component of spin in part A, as 
indicated in the inset. Eor OES, the subsystem A contains /a = 4 orbitals, where we show data for Na = 6 and Na = 7. The 
spectrum is resolved as a function of momentum of subsystem A. 


I = 3, n = 2. We have suppressed, as usual, the spinor 
and Gaussian parts of the wave function. The permanent 
was defined in Eq. (42). 

We are particularly interested in the states 
and ^ 221 - The former was introduced by Moore and 
ReacP, and subsequently analyzed in detail by Read and 
RezayP^ (see also Ref. EOl). As a prototypical example 
of a state that derives from a non-unitary GET, it was 
shown to describe a critical point between the ferromag- 
net and paramagnet. Its thin-torus limit was recently 
studied in Ref. |39l The latter state, ^ 221 , has been ad¬ 
dressed by Ardonne and RegnaulfP^, who computed some 
of its properties by identifying its root configuration on 
the sphere and deriving its “squeezing” properties. Here 
we show how to write down the parent Hamiltonian for 
these two states for the cylinder and torus geometry. This 
is of particular interest in the case of ^ 221 , where such 
a Hamiltonian has not been reported for any geometry 
before. 

To begin with, we want to derive many-body projec¬ 
tion Hamiltonians whose null space contains We 

first consider 2-body interactions. From Eq.[^ van¬ 
ishes as the power when Zi Zj (or Wi Wj). Hence 
we need 2-body terms where m < I and |5'| = 1, 

i.e. in the |tt) oi* \ii) channels. In the S' = 0 (or |ti)) 
channel, however, one term in the permanent will cause 
^iin not to vanish for Zi Wk- Hence we do not need 
any C^^,5=o- 

For 3-body interactions, we find that the Halperin 
function parts contribute a degree of / + 2n for S = ^ 
(or the 6tc.) channels. Also, the perma¬ 

nent will remove one degree. Hence we need 
m < / -f 2 n — 1. An analogous analysis for |ttt) oi* liii) 
reveals that we also need 3^2 m <31. The matter 

simplifies a bit since not all of the requisite U'^ 


2,3 

m,S 


men¬ 


tioned above actually exist. In Sec. IV and also Table 1 


of Ref. [43l the PPs g are nonzero for even m when I 
is even and S = 1, or / odd and S = 0 (and vice versa 
for total anti-symmetry in spin-orbit space). For A = 3, 
1/2 6^ists for m > 1, but 3/2 exist for 

m = 0,1, 2 in the case of fermions (odd /). In summary, 
the above considerations for n = 1 and / = 1,2 imply 
that only is required to produce ^m, consistent 

with Ref. [TTl On the other hand, Vq 1, ^21/2’ 

and U^^j 2 for m = 0,2, 3,4, 5 are all required to pro¬ 
duce ^221- Therefore, we see that formulating the par¬ 
ent Hamiltonian for ^221 is indeed rather involved, and 
crucially benefits from the systematic approach we follow 
here. 

For completeness, we quote the explicit expressions 
for these two Hamiltonians in the form = 

+ (It) ^ It))- For the state, we have 

(56) 

n 

with the corresponding operator and polynomial ampli¬ 
tude 

6y/2 = K^/2 ^ Pi.i/2(K;{n})e“'^^"c„iC„2C„3, 

ZI ni=n 

(57) 

Pi.1/2 = («i - ^2) Itti) + (^2 - na) 14-tt) 

+ (fia - fii) mt), (58) 

where as before hj = n, — | Yhj nj- For $ 221 , we have 
H 221 = Uli + Z 7mf7^,i/2 

m=l,2 

T ^ 7 n^m, 3/25 

m= 0 , 2 , 3 , 4,5 
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FIG. 10. (Color online) PES (left) and OES (right) for the ^221 state on a finite cylinder with aspect ratio equal to 1. Eor PES, 
we consider iV = 8 bosons and the subsystem A contains Na = 4 particles. The PES is plotted as a function of momentum 
Ky of subsystem A. Different values of S^, the total 2 ; component of spin in part A, are indicated in the inset. Eor OES, we 
consider iV = 10 bosons, the subsystem A contains Ia — ^ orbitals and the OES is plotted as a function of momentum Ky of 
subsystem A, for two values of the total number of particles in subsystem A (Na = 5 and Na = 6). 


where again 7 ^ and r]m can be chosen to be arbitrary 
positive constants. 

We have numerically diagonalized the Hamiltonian 
(59) and verified that it gives a unique zero energy ground 
state on finite cylinders when = 3N/2 — 3. By stretch¬ 
ing the cylinder, we confirmed that the ground state re¬ 
duces to the pattern (t, i) 00 (t, i) 00 ..., which coincides 
with the root configuration given in Ref. |64] In Fig. 
(left) we compute the PES for N = S bosons, which can 
be directly compared with Fig. 7 in Ref. |64| As ex¬ 
pected, we find an agreement in the counting of the PES 
levels between our results on the cylinder and the sphere 
data in Ref. [64l Finally, we also provide a plot of the 
OES in Fig. (right). As shown by Li and Haldan^^, 
a non-Abelian state like the Moore-Read state will have 
different counting depending on the position of the or¬ 
bital entanglement cut, which effectively probes different 
topological sectors of the theory. This is also the case 
in our example, as we see the counting changes depend¬ 
ing on the number of particles in subsystem A while the 
location of the cut remains fixed. Furthermore, the spec¬ 
trum is fully chiral, which suggests that this state may 
be described by a simple CFT as conjectured in Ref. [64l 
The CFT can in principle be identified using the counting 
of the OES in Fig. 10 For this purpose, larger systems 
may be required in order to verify that the counting we 
see in Fig. is saturated and indeed corresponds to the 
thermodynamic limit, which we defer to future work. 


VI. CONCLUSION 

We have provided a general method for constructing 
quantum Hall parent Hamiltonians with the desired type 
of clustering properties on cylinder and torus geometries. 
The approach developed here completes the program ini¬ 
tiated in Ref. uni which is suited for the disk and sphere 


geometry that have a conserved ^-component of angular 
momentum. We have performed extensive analytic and 
numerical checks of the proposed model Hamiltonians, 
and generally find complete agreement with previous re¬ 
sults if available. We have also demonstrated that it is 
possible to construct the parent Hamiltonians for some 
rather complex spinful non-Abelian states which previ¬ 
ously have not been provided in the literature. 

The method presented here opens up several future di¬ 
rections. One can use the Hamiltonians derived here to 
systematically scan for filling factors where the Hamil¬ 
tonians have zero-energy ground states, in the spirit of 
Ref. |671 Performing this kind of search is much more nat¬ 
ural in the torus geometry which does not suffer from the 
known shift bias^. Furthermore, our approach straight¬ 
forwardly generalizes to 4-body interactions whose possi¬ 
ble incompressible ground states have not been systemat¬ 
ically investigated before. Another interesting extension 
to pursue is the study of multicomponent states beyond 
the familiar SU( 2 ) spin, such as the SU(4) valley/spin 
symmetric non-Abelian states. They may have relevance 
in the realization of FQHE in graphene, where the open 
accessibility of the electron gas may a llow for further tun- 
ability of the interaction profild^^^. 


ACKNOWLEDGMENTS 

CH Lee is supported by an NSS scholarship from the 
Agency of Science, Technology and Research of Sin¬ 
gapore. ZP acknowledges support by DOE grant de- 
sc0002140, Research at Perimeter Institute is supported 
by the Government of Canada through Industry Canada 
and by the Province of Ontario through the Ministry of 
Economic Development & Innovation. RT is supported 
by the European Research Council through ERC-StG- 
TOPOLECTRICS-336012. 








19 


Appendix A: Barycentric coordinates for many-body 
pseudopotentials 

We describe how to find a manifestly S'^r-symmetric 
embedding of the tuple (ni,riAr) onto the {N — 1)- 
dimensional simplex in This can be achieved in 

barycentric coordinates, which is also useful for diverse 
applications involving permutation symmetry with a lin¬ 
ear constraint!^. The most straightforward construction 
is to write 


N 




(Al) 


k=l 


where x G and {/3k}, k = A forms a linearly 

dependent set of basis vectors, also in normalized 

so that 

- - A 1 

f3j ' ~ 


N-1 


N -1 


(A 2 ) 


Geometrically, the define the vertices of a simplex, 
and are at angle of cos“^(— 1 /(A — 1 )) from one another. 
The vertex k corresponds to the least isotropic configu¬ 
ration with fik oc and fij (x ^ j ^ k. A basis 

consistent with the above requirements is 


/3i = (1,0,...,0), 


(A3a) 

/32 = (C'i,5i,...,0), 


(A3b) 

/33 = (Ci,5iC2,5i52,...,0), 


(A3c) 

04 = {Cl, 8,02,8,8203,818283,. 

.., 0 ), 

(A3d) 

,., = {C,,8,C2,8,82C3,...,8,.. 

• Sn-2), 

(A3e) 

0N = {01,8,02,8,8203,...,-8,. 

' • • Sn-2), 

(A3f) 


with + = 1, 1 < k < N — 2. Upon enforcing 

the scalar product constraint in Eq. |A 2 [ we require that 
Sk+I = 1 - (1 + l/^)/(nj=i which also implies 

that (7| + S‘f.Ck +1 = Ck‘ A simple solution fortunately 
exists: 

c. = (A 4 ) 

which implies that the k-th. projected component of the 
relative angles between the position vectors of the ver¬ 
tices approaches tt as /c and N increase. We proceed 
by substituting Eq. |A4| into the explicit form of vertex 
positions l3k, and expressing the latter in terms of the 
spherical coordinates. Erom Eq. |A 1 [ we can easily check 

N-1 


rik = 


Nk 


-X’ Pk, 


(A5) 


and that 

= Ifp 


Nk 

N~k 


i5:( 


n y 
Ay 






(A6) 


Eor the purpose of orthogonalizing the PPs over the inner 
product measure in Eq. we will also need to express x 
explicitly in terms of angles in origin-centered spherical 
coordinates: 


xi = W cos (/:?!, 

(A7a) 

X 2 = W sin ifi cos {p2 5 

(A7b) 

N-3 


xn- 2 = w cos cpk sin(/ 9 /c. 

(A7c) 

k=l 


N-2 


xn -1 = W sinipk. 

(A7d) 

k=l 



Substituting the explicit expressions from Eqs. |A3| and 
|A7| into Eq. |A5[ we obtain 


(A 8 ) 


w 


n 2 = cos (fii + ^N{N -2) sin (pi cos <^ 52 ) , 

(A9) 


N 


sm ipi cos ip 2 


y I m (nLiSin</?i) cos(^j+i 

■EV(iV-j-2)! {N-j){N-j-l) 
N\ (nti^ sin cos (pk 


j=2 


+ ■ 


{N-k-l)\ N-k + 1 


(AlO) 


W 


= V7V-2 


N 


sm (pi cos If 2 


I m (riLi sin V’i) cos fj+i 
{N-j-2)\ (N-j){N-j-l) 
rN-2 \ 

(All) 





sm Lpi cos (^2 


I TV! (OLi sin ‘Pi) cos fj+i 
{N-j-2)\ {N-j){N-j-l) 
m-2 \ 

(A12) 



These are the explicit expressions for transcribing the 
fij^ 1 < j < N indices directly into spherical coordinates. 
In (AlO), k ranges from 3 to A — 2. 
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Appendix B: Second-quantized pseudopotentials vs. 
real-space projection Hamiltonians 

We provide a brief comparison between the second- 
quantized PPs that have been derived in this paper and 
real-space projection Hamiltonians elsewhere in the lit¬ 
erature, e.g. Trugman-Kivelson type Hamiltonians on 
the infinite plane. We show that our geometric PP con¬ 
struction avoids certain ambiguities that plague the latter 
approaches. 

Neglecting internal DOFs for simplicity, a generic real- 
space PP living in total relative angular momentum sec¬ 
tors up to m takes the form 

j - f2)...S^{fN-l -Tn) 

(Bl) 

such that ^jdj = m. The various dj’s refer to how 
the derivatives are distributed among the particles. In 
this form, there is no simple one-to-one correspondence 
between the dj values and the relative weight of H in 
the various m sectors. To find the relative weights, one 
has to project H onto the sector spanned by a state with 
angular momentum m. It is most convenient to adopt 
the symmetric gauge, where such a state takes the form 

l^'m) OC \ zi ,..., zn ) where p is a 

symmetric or antisymmetric polynomial and Zj = Xj-\-iyj 
refers to the particle position {z defined below): 


where = ^dzdz has been used. The delta func¬ 
tions, which enforce that all Zj be equal, also enforce 
^3 ~ ^3 ~ ^ ^3 — j' Thus we only get a nonzero 

contribution from constant integrands. Specific examples 
were worked out in the appendices of Ref. [42l in general, 
an operator H has an overlap with various m sectors, and 
a PP that lies purely in one m sector must be a compli¬ 
cated linear combination of the Trugman-Kivelson terms 
with various sets of dj’s. Some progress was made in 

Ref.[44l where operators containing Lm (^2 (iv^l)vf ) 
shown to be contained purely in one m sector. That, 
however, does not always work for fermions as will be 
explained below. 

The mathematical complications above are further ag¬ 
gravated when the space of PPs in the sector m is de¬ 
generate. In particular, the matrix is not 

of full rank and some combinations of bona-fide PPs of 
degree m will still not give a positive energy penalty to 
states in the m sector^. In particular, it is possible for 
certain Trugman-Kivelson Hamiltonians with 2m deriva¬ 
tives to evaluate to zero upon taking into account the 
fermionic antisymmetry. Consider for instance the N = 3 
fermionic state IT 3 ) oc ( 5 i- 52 )(^ 2 -^ 3 )(^ 3 -^i)ki, ^ 2 , ^s)- 
The following operators have a zero projection on it: 

(^'3|V®(52(ri - r2)S\r2 - r3)|«'3) = 0 

and 


(n/.v, 


(B2) 


Upon substituting and integrating by parts, we can 
express as 

j “ f2)...<5^(fAr_i - rjv) 

X (B3) 


which can be evaluated as 


(^'3|V?ViVi52(ri - r2)5\r2 - r3)|^'3) = 0. 


Indeed, the only nonvanishing Trugman-Kivelson type 
fermionic PP in the sector is — 

r 2 )d‘^{r 2 — rs). As jTs) is of degree < 2 in each variable, 
H must contain a derivative of degrees 2dj < 4 in each 
variable for a nonzero projection. Since only constant 
terms in the integrand of Eq. |B4| contribute, we also 
need derivatives in all three particle coordinates. 

By contrast, our PP construction appeals to the or¬ 
thogonality structure of the PPs from the outset^ and 
avoids all of the above problems that generically arise 
when one attempts to use the Trugman-Kivelson ap¬ 
proach in more complicated many-body cases. 


Appendix C: Pseudopotential parent Hamiltonian 
from a given wave function 


: 4' 


^ ^IJ J dzjdz*^ 5'^{zi - Z2 )...S^(zn-i - zm) 

U=i ' / Vi=i / 

(B4) 


We briefly summarize the heuristic procedure how to 
determine the PPs whose null spaces contain a given 
QH state. These are the PPs that should be included 
in its parent Hamiltonian, since they penalize denser or 
“simpler” states that will otherwise be realized at the 
same filling fraction. The method consists of two steps: 
(i) identify the thin-torus root pattern from the first- 
quantized wave function (Sec. Cl), and (ii) use th e roo t 
pattern and power counting to obtain the PPs (Sec. C 2). 





21 


1. Wave function to root pattern 

Suppose we are given a QH wave function with the 
polynomial part ^2, ..O? where the z^’s denote the 

positions of particles. In the infinite plane, Zj is the usual 
complex coordinate Zj = Xj + iyj. Since we will work 
on a cylinder geometry, however, we need to perform a 
stereographic mapping: Zj , where n = 211 jLy is 

shown in Fig. For simplicity, assume there are no 
internal degrees of freedom. 



-<- 

H = 2'KlBlLy 


FIG. 11. (Color online) In the limit of thin cylinder {n 
oo), the single-particle orbitals become well-separated from 
one another. As a consequence, the dominant configurations 
(“root patterns”) in the expansion of a many-body state ^ 
are those Fock states where particles minimize their classical 
electrostatic energy. 

To each wave function defined on a cylinder, we can as¬ 
sign one or several thin-cylinder root patterns. Those are 
the Fock states with non-zero weight as the circumference 
of the cylinder is taken to zero {k oc). Typically, any 
given FQH state \^p) in the thermodynamic limit decom¬ 
poses onto many Fock states, as \ 2 p) represents a strongly 
correlated state. However, when the cylinder is stretched, 
the weights of most of the Fock states in the decomposi¬ 
tion vanish. The only Fock states surviving in the limit 
hz ^ oc are called “root patterns” and are useful for 
finding the parent Hamiltonian (Laplacian) for as its 
kernel. 

In order to find the root pattern, we must examine the 
dominant terms in the decomposition of t/; as /i: ^ oo. 
Usually, this amounts to solving a classical electrostatic 
problem. Here, we consider the case of an open cylinder, 
where the root pattern will typically be unique. To get 
the remaining root patterns, that arise due to topological 
degeneracy or non-Abelian statistics, one must consider 
the system on a torus. Due to the complicated form of 
torus wave functions, it is more convenient to work on a 
cylinder but perform a flux insertion which allows one to 
access other topological sectors. Details of this are given 
e.g. in Ref. [56| 

Fixing the total number of particles, we study the 
dominating term in It is z^^z'^'^z^^... such that no 
other term (denoted by ') has m[ > mi or, > m 2 
if m'l = mi or, > m3 if m'l = mi and m '2 = m2, 
etc. The root pattern can then be written down as the 
partition defined by mi, m2, 

(mi, m2, m3,...) ^ 10...010...01... (Cl) 

where the Us are at position mi,m2,.... If two m^’s 
coincide, as can occur for bosons or multicomponent 


fermions, we shall label that position by ’2’ and so forth. 
For instance, consider the 1/m Laughlin state with 2 par¬ 
ticles. We have i^{zi^Z 2 ) = {zi — Z 2 )^ which, when ex¬ 
panded out, contains monomials 

(C2) 

Evidently, the term z^^z^ is in the root pattern. Had 
we considered a droplet of N particles, this term would 
have become with total de¬ 
gree and the root 

(C3) 


where 0"^“^ represents a string of m — 1 consecutive ze¬ 
roes. For this simplest case, a droplet of N = 2 particles 
is sufficient to determine the root configuration. Because 
we deal with translationally invariant states, a root con¬ 
figuration for larger N is obtained by repeating the min¬ 
imal pattern. 

Let us now consider a more complicated example such 
as the generalized Pfafhan state introduced in Eq. (46) 
in Sec. By using the identity [Pf{A)]‘^ = Det(A^^ 
one can show by induction starting from n = 4, the 
most compact particle droplet, that the dominating term 
is z^^~^Zo^z^~^. This defines the root configuration in 


Eq. (47), 10^“^ 10^ 10^“^ 10^.... It is also easy to read off 
the filling fraction from the above: we have 2 particles 
per period, and each period has length 2 q. Hence the 
filling fraction is ^ = 1/g. 


2. Root pattern to parent Hamiltonian 

Given a root pattern, we can easily infer what rela¬ 
tive angular momentum states should be penalized to 
construct a projector onto the state. Eor instance, if 
the closest pair of particles are ...10^1..., no two particles 
have relative angular momentum less than / + 1. Hence 
the state should lie in the null space of any two-body PPs 
with m < 1. Consequently, the parent Hamiltonian 
should be a linear combination of these PPs, so that it 
penalizes all states denser than the given one. This is 
discussed at length in Ref. [30l 

To figure out the required many-body PPs, we apply 
the above procedure analogously to clusters of particles 
with minimal relative angular momentum. Note that 
some many-body PPs need not be included if the inter¬ 
action is already precluded by particle (anti)symmetry: 
Consider the 1/q Pfafhan state described above, with a 
root given by Eq. Obviously, no two particles have 
relative angular momentum less than q — 1. Naively, we 
will then need to include all with m < q — 2. The 
1/q Pfafhan state, however, is fermionic (bosonic) when 
q is even(odd), and Ug_2 exactly vanishes in both cases. 
Hence the parent Hamiltonian only contains U^<g_2- 

Let us hnally analyze 3-body interactions. Erom the 
root pattern, the minimal total relative angular momen¬ 
tum is / = {q — 1) 2q = ‘iq — 1. However, conhgurations 
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with / = 3g — 2 = 2(g — l)H-g automatically vanish since 
they require pairs of particles with relative angular mo¬ 
menta q and q — 1 simultaneously, which are forbidden by 
either fermionic or bosonic statistics. With lower 3-body 
relative angular momentum already excluded by the 2- 


body PPs, the only 3-body PP we need to include in 
the parent Hamiltonian is pattern re¬ 

peats from the third particle, there is no need to consider 
interactions with 4 or more bodies. 
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